Two photon decay of Z' as a probe of Bose symmetry violation at the CERN LHC 
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The question if the Bose statistics is broken at the TeV scale is discussed. The decay of a new 
heavy spin 1 gauge boson Z' into two photons, — >■ 77, is forbidden by the Bose statistics among 
other general principles of quantum field theory (Landau- Yang theorem). We point out that the 
search for this decay can be effectively used to probe the Bose symmetry violation at the CERN 
LHC. 
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I. INTRODUCTION 

For all known particles, there is a remarkable one-to- 
one correspondence between their spin type and statistics 
type: 



integer spin o commutation relations 
exchange — symmetric wave functions 

Bose statistics; 



(1) 



half — integer spin ■(-> anticommutation relations o 

exchange — antisymmetric wave functions -H- (2) 

Fermi statistics. 

In quantum field theory, this spin-statistics connec- 
tion (SSC) can be understood in different ways. For free 
fields, the connection was established using conventional 
lagrangian and group-theoretical methods ^1] . 

Interacting fields were tackled using axiomatic quan- 
tum field theory (without any lagrangians) , and by early 
1960s the celebrated spin-statistics theorem was proved 

It is worth recalling that the axiomatic treatment of 
the spin-statistics connection Q does not cover the case 
of quantum electrodynamics and other gauge theories in- 
cluding the Standard Model (for a detailed explanation, 
see (3, Sec. 8.1). 

The physical root of this difficulty is the fact that while 
the photon is a massless spin-1 particle with only two 
polarizations, the corresponding 4-potential yl^ has Jour 
components. 

Consequently, ordinary Hilbert space with positive 
metric cannot accomodate the 4-component field oper- 
ator A^, and the introduction of indefinite-metric space 
is required. 

This is at odds with the axiom underlying the proof of 
the spin-statistics theorem, which states that the Hilbert 
space metric must be positive. 

As a result, the original set of axioms Q needs to be 
modified to include Hilbert spaces with indefinite metric 
as well. Although progress has been achieved in this area 
(see, e.g., 0], Ch.lO) the extension of the spin-statistics 
theorem to gauge theories is yet to be formulated. 



In any case, the theorem does not forbid small vio- 
lation of the SSC and the question whether it exists or 
not remains open. (That does not mean, however, that 
construction of a theory with small statistics violation is 
easy. To date, most attempts to find a local relativistic 
quantum field theory with small statistics violation have 
been unsuccessful.) 

There have been several works that have attempted 
to improve the original spin-statistics theorem by going 
beyond the Bose-Fermi alternative 4], and, in particular, 
to rule out the small violation of SSC. However, they also 
assume positivity of Hilbert space metric and hence do 
not apply to gauge theories. 

Even if the small SSC violation was shown to be the- 
oretically forbidden on the basis of general principles of 
QFT such as Lorentz invariance, locality etc., the experi- 
mental tests of SSC would still be of interest because they 
would be important tests of those general principles. 

We may assume that at some extremal or small enough 
distances where the usual notions of the local fields be- 
come invalid, their proper reconsideration and a general- 
isation should be required, which can bring us the possi- 
bility of changing or correcting the usual Spin-Statistics 
Correspondence. 

Since 1987 there have been significant theoretical and 
experimental efforts to motivate and find tiny departure 
from the established connection between spin and statis- 
tics. 

Originally, most efforts [5-30], especially in the exper- 
imental field, were actually concentrated on discussing 
small violation of the Pauli exclusion principle rather 
than violation of Bose statistics. Many dedicated ex- 
periments have been performed to give strong bounds 
on the violation of the Pauli principle. Also, the topic 
was discussed in the context of string theories [3l| and 
cosmology [32,]. For recent reviews, see [3^. 

Later, Bose statistics came under scrutiny as well. Ini- 
tially, experimental proposals of searching for deviations 
from Bose symmetry used the spin-zero nucleus of oxygen 
1^0 as the test object [Mill. 

For a review of subsequent experiments with O2 and 
CO2 molecules, see [36]. 

The first experimental upper limit on the validity of 
Bose statistics for photons was obtained in Ref. [13] 
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based on the idea of forbidden two-photon decays of Z- 
boson. 

The same idea was later explored in Ref. 38], in the 
context of low-energy, high-sensitivity atomic two-photon 
transitions (see also |39l|). 

This approach exploits one of the important conse- 
quences of the Bose symmetry, first observed by Lan- 
dau ^] and Yang [41]: a pair of photons cannot form a 
state with total angular momentum equal to unity. The 
Landau- Yang theorem uses general principles of rotation 
invariance, gauge invariance, and Bose statistics to de- 
rive certain selection rules for decays of a parent particle 
into two photons. For a parent with spin one, the de- 
cay amplitude into the exchange-symmetric state of two 
photons vanishes. 

Therefore, the decay of any spin-1 boson into two pho- 
tons is absolutely forbidden (for a textbook proof, see, 
e.g., 42]). If, however photons do not obey Bose statis- 
tics, there will be a nonzero decay amplitude involving 
two photons in an exchange-antisymmetric state (37| . 
This provides a clear reason why the diphoton system 
is especially interesting in testing the degree with which 
Bose symmetry is exact. 

Recent experiments have explored the possibility of 
small violations of the usual relationship between spin 
and statistics, which are impossible within conventional 
quantum field theory. Experiments searching for the 
transitions between atomic states with J = and J' = 1 
for degenerate photons (i.e. photons of equal energies) 
test the Bose statistics at the eV scale and yield up- 
per limit on the ratio v of the rate of statistics- violating 
transitions to an equivalent statistics-allowed transition 
rate, u < A.Q x 10"" at the 90% confidence level p3l| . 
At higher energies the limits on the branching frac- 
tion of two-photon decays of the triplet positronium (or- 
thopositronium) Br{oPs 77) < 2.4 x 10~^ [13] and 
charmonium Br{xci 77) ^ 3.5 x 10~^ jl^ have been 
reported. At the energy scale of 100 GeV the results of 
the search for the gauge boson Z into two photons were 
obtained at LEP ^ Br{Z 77) < 1.4 x IQ-*. 

From general arguments, one can expect that any vio- 
lations of Bose symmetry, if any, would be better mani- 
fested at higher energy scales. Then the question arises: 
could interesting physics be found by combining the idea 
of searching for Bose statistics violation and assumption 
that Bose symmetry might be broken at a high energy 
frontier at the CERN LHC? Out of aU (neutral) spin 1 
bosons it is natural to concentrate on the heaviest one- 
the new heavy gauge boson Z' which appears naturally 
in many extensions of the SM. 

The explicit goal of this paper is to extend the results 
of the previos work [37] to the TeV scale and to show 
that the search for the decay Z' — 77 at the CERN 
LHC could result in a more radical departure from stan- 
dard physics: the possible observation of a small violation 
of Bose statistics, which would provide strong evidence 
for the existence of new physics. The rest of the paper 
is organized as follows: in Sect. II we describe the phe- 



nomenological model (37[ of the Bose symmetry violation 
and its extensions; we then show how the parameters of 
the model are constrained by theoretical arguments and 
clectroweak precision data. Following this. Sec. Ill deals 
with the Z' sector in some detail. The main characteris- 
tics of the Z' boson are described within the framework 
of several Z' models. Section IV- VI present the results of 
a sensitivity study for the Bose symmetry violating pro- 
cess pp^ Z'^ 77 at the CERN LHC at = 14 TeV. 
Finally, Sec. VII presents the conclusions of this work. 

II. A MODEL FOR BOSE SYMMETRY 
VIOLATION 

In this section we describe a simple model of Bose sym- 
metry violation suggested in Ref. |37| and its extensions. 

The method '37'] is to write down the most general 
form of the decay amplitude of the spin-1 particle into 
two photons and then apply the conditions of gauge in- 
variance and Bose symmetry to that amplitude. If both 
conditions are applied, the resulting amplitude is exactly 
zero. However, if we impose the condition of gauge in- 
variance but do not require the Bose symmetry, the re- 
sulting amplitude is not zero. We then can obtain the 
two-gamma decay rate of Z-boson and compare it to the 
experimentally known upper bound on the branching ra- 
tio oi Z ^ 77. In this way a direct bound on Bose 
symmetry violation for photons can be obtained [s?]. 

Now, a few remarks about the relation of this method 
to alternative models of small Bose symmetry violation. 
The most well-studied model is "the quon model" pro- 
posed in [13]. Quons are particles described by the com- 
mutation relations of the form: 

aka'^ - qa^lak = Ski- (3) 

However, it was shown in Ref. [48] that in relativistic field 
theories quons must be either fermions or bosons. For 
this reason, the quon theory cannot be used for the anal- 
ysis of the decay Z — > 77. 

Let us turn now to the construction of Z — 77 decay 
amplitude. We require that this amplitude satisfies all 
the standard conditions, such as relativistic invariance 
and gauge invariance, but we do not require this am- 
plitude to be symmetric under the exchange of photon 
ends. The most general Lorentz invariant form of the 
amplitude S of the decay Z — > 77 is: 

5(fci,fc2,ei,e2) = CA^z.(fci, fc2)eoeie2 (4) 

where fci and /c2 are photon momenta, ei and £2 are 
photon polarization vectors, eg is Z-boson polarization 
vector. 

Even though violation of SSC could require viola- 
tion of Lorentz invariance as well, using more general 
parametrization in Eq. @ would be overkill. 

Note that terms in ca/^i/ proportional to fcip,/c2i/ and 
^lA + do not contribute to S due to the conditions 

e^fci^ = 0, e^fe. = 0, e^(fciA + k2x) = (5) 
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and can therefore be ignored. 

We focus first on the part of ca^i/ that does not contain 
the e-tensor: it has the following Lorentz invariant form: 

ca^i. = bi9Xf^ki,y (6) 
+b2g\^k2f, + gg^_iu{kix - k2x) + h{kix - k2\)ki^k2f. 

Now, the condition of the electromagnetic gauge in- 
variance reads 

c\fj.uki = 0, c\^vk^ = 0. (7) 
These conditions can be satisfied if we put 



6i = 62 = 0, h = - 



Then the amplitude cx^y becomes [37| 



cx^v = g{ki- k2)x{gi_iu 



Ml 



(8) 



(9) 



In principle, g could depend on some scalar products 
of the momenta, but in our case, since all the particles 
are on mass shell, we have fci/c2 = Ml/2 (and, of course, 
k\ = kl = G), so that g is a pure number. Note that 
the above amplitude automatically satisfies the condi- 
tion (fci + k2)^c\i^i^ — 0. We see that this amplitude, as 
expected, violates Bose symmetry because 



CAAiiy(fcl, ^2) = -CX,^fj.ik2,ki) 



(10) 



whereas Bose symmetry requires cxp,^{ki,k2) — 

+CAi/Ai(fc2, fci)- 

Thus the parameter g can be interpreted as the param- 
eter of Bose statistics violation which will be marked by 
the subscript Bl: g = qbi- 

Next, it can be shown that the following Lorentz- 
invariant terms containing the e-tensor also satisfy the 
conditions of gauge invariance and Bose-antisymmetry: 

cx^i.v = UlC'^i^ ^ k2x)efj,^ai3k"k§ (11) 



and 



CXtJ.u = ^^[{^Xf^afskli, - exvai3k2,i)k"k2 + (12) 



(-XiJ.va{kl + fc2)"(fcifc2)]. 



Now, calculating the width of the decay Z — >■ 77 with 
the help of the amplitude Eq. (|9]) we obtain [t^ : 

ri(^^77) = T7^|5P = 4^. (13) 



247r 



Exp erimentally, it has recently been measured at LEP 
[49| that 



Br{Z 77) < 5.2 X 10" 



(14) 



Therefore 

ri(z^77) 



gliMz 



< 5.2 X 10' 



(15) 



Thus, taking into account that TtotiZ) ~ 2.5 GeV !49|, 
finally, we can obtain our upper bound on the Bose vio- 
lating coupling 



gsi < 1-1 X 10 



(16) 



We see that in the framework of this model the rate of 
the decay Z — >■ 77 is small; however, it can be enhanced 
for higher-mass particles. Hence, the heavier spin-1 Z' is 
a good candidate for the searching for effect of Bose sym- 
metry violation through the Z' decay into two gammas. 

A similar analysis can be carried out for the ampli- 
tudes of Eq.lfTTjl and (fT2l) . The decay widths due to the 
amplitudes (|llll2p are equal to: 



r2{Z 77) 



967r 



. r3(Z^77) = ^. (17) 
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The corresponding upper bounds on the constants 
gB2, gB3 are: 



gB2, gB3 < 2.2 X 10' 



(18) 



In the context of our approach, the reason for treating 
the three amplitudes separately is purely technical and 
does not involve the essential physics. Indeed, adding the 
three amplitudes with three arbitrary coupling constants 
would be possible but it would only add considerable 
complexity and obscurity into experimental simulations 
without any gains in physical understanding. 

A few remarks are now in order concerning the rela- 
tionship between the approaches developed in Ref. [s^l 
and a later work, Ref. ^38,] . 

In Ref. [11] it was claimed that the limit obtained in 
[S^l is too weak to be of any significance. This conclu- 
sion was reached within a very specific model for Bose 
symmetry violation which is different from the model- 
independent approach suggested in [s^l- A detailed dis- 
cussion of similarities and differences between the two 
frameworks would be out of place in the present paper, 
so we restrict ourselves to a few general comments only. 

In a nutshell, the difference is this: Instead of con- 
structing the most general Bose-violating amplitude of 
the decay and the coupling constant 175, the authors of 
Ref. [3^ introduce the "probability for two photons to be 
in an antisymmetric state" v. The physical decay width F 
is then obtained as F = Tb + i'Tp, where F^ is the width 
for the decay into ordinary "bose-photons" , and Fj? is 
the width of Bose-symmetry violating decay into "fermi- 
photons". For instance, if the decay amplitude for the 
usual, "bose-photons" 1 and 2 is = A12 + ^21, then 
the corresponding amplitude for "fermi-photons" would 
be Ap = A12 - A21, with Tb oc \Ab\'^ and Ff oc \Af\'^. 

This approach may look simple and natural, but five 
points should be kept in mind: 
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1. Well-known difSculties arise when v is introduced. 
According to quantum mechanics, every probability is 
the square of (modulus of ) amplitude. So, after we in- 
troduced V we must introduce the "amplitude for two 
photons to be in an antisymmetric state" , let us call it 
[1. 

So, the Bose-violating two-photon state will be \B) = 
iVdS*) -I- {N is the normalization factor). Next, 

due to superposition principle, we can add to it the state 
-A^IS'), obtaining \B') = Now, \B') is obviously 

not properly normalized, and after normalization, it be- 
comes just \B") = i.e. "amplitude for two photons 
to be in an antisymmetric state" becomes 1 instead of [i. 

So it is hard to give physical meaning to the concept of 
"probability for two photons to be in an antisymmetric 
state" , unless superposition principle is modified in some 
way. 

2. It is believed [11] that the 'Bose' and 'Fermi' ampli- 
tudes Ab and Ap not interfere, i.e., the total proba- 
bility is assumed to be F oc \Ab\^ + v\Ap\^ rather than 
F (X I As -I- [iAp\^ . In Ref. [38] this assumption is jus- 
tified by invoking the rule that the matrix elements of 
a symmetric Hamiltonian has zero matrix elements be- 
tween the states of different symmetry Q. However, this 
rule appears to be superfiuous as the model [111 itself 
dictates whether Ab and Ap interfere or not. Indeed, 
in this model, the necessary and sufficient condition of 
non-interference is \A\2\ = ]^2i] (assuming for simplicity 
that jjL is real). (It should be noted that whether the su- 
perposition is coherent or incoherent is irrelevant for the 
decay Z' — 77, which is forbidden at zeroth order.) 

3. Should we consider 1/ is a 'universal' parameter, i.e., 
independent of the physical process, energy scale etc.? 
This may or may not be true depending on the details 
of the underlying specific theory. For instance, it is not 
inconceivable that v could be energy-dependent and grow 
with energy. In this case, upper limits on v obtained 
in low- and high-energy processes would not be directly 
comparable to each other. 

4. Finally, Ref.[3l] makes (implicitly) a strong but ar- 
bitrary conjecture that the amplitude A12 should be cal- 
culated in the Standard Model. There is nothing wrong 
with it as long as we remember that: 

— this is just one of many possible assumptions, and cer- 
tainly not the most general or unique. 
— the parameter v is meaningful only if this assumption 
is made. 

5. If we pursued a similar approach, the result would 
be that the rate of Z' — >■ 77 decay becomes rather small: 
due to a fermion loop it would be suppressed by at least 
a factor of (cf. [38]). As a result, its search at LHC 
would require a higher sensitivity. 

By contrast, Ref. [13] did not rely on this or similar 
conjectures, but tried to be more general. In this general 
approach, there is no point in estimating v, because this 
parameter belongs to a different model based on addi- 
tional specific assumptions. 

We feel that at the moment our level of understanding 



is, unfortunately, insuflicient for telling on purely theo- 
retical grounds which approach is correct, and we need 
an experimental search that could settle the issue. 



III. THE Z' DECAYS 

Consider now the new boson Z', which appears in 
many models of physics beyond the SM, see e.g. [50j-[58|. 
The Z' is assumed to be a more massive than the gauge 
boson Z of the standard model. The most direct chan- 
nel to probe the existence of the Z' at a hadron collider, 
such as the CERN LHC [111, is the Drell-Yan process. 
The Z' that decay to leptons, pp ^ Z' ^ -\- X with 
I = e, fi, have a simple, clean experimental signature, and 
potentially could be discovered at the LHC with a mass 
up to 5 TeV, see e.g. [13, HH. This new object is sup- 
posed to be neutral, colorless and self-adjoint, i.e., it is 
its own antiparticle. The mass of the new boson could 
be identified unambiguously by a study of a resonance 
peak in the dilepton invariant mass distribution. The Z' 
may be classified according to its spin, which could be 
defined by measuring the dilepton angular distribution 
in the reconstructed Z' rest frame. The Z' could be a 
spin-0 V in R-parity violating SUSY, a spin-2 Kaluza- 
Klein (KK) excitation of the graviton as in the Randall- 
Sundrum (RS) model, or even a spin-1 KK excitation of 
a SM gauge boson from some extra dimensional model 
[53| . Another possibility for the spin-1 case is that the Z' 
is a true Z', i.e. a new neutral gauge boson, which is the 
carrier of a new force, arising from an extension of the 
SM gauge group. For much more extensive discussions 
of specific models and other implications see several ex- 
cellent reviews (521-I53. [stI |62| , and a more complete list 
of references therein. 

The current best direct experimental lower limits on 
the mass of Z' of a few popular models came from the 
Tevatron and restrict the Z' mass to be greater than 
about 900 GeV when its couplings^ to SM fermions are 
identical to those of the Z boson [49| . 

Consider now the allowed branching fraction Br{Z' — > 
77) — r(Z' — > 77)/F(Z' — > in several interesting 

Z' models. First, we shortly describe these models and 
the SM fermions couplings the Z' . 

• the Eg models are described by the breaking chain 
Eq SO{l{))y.U{l)^ SU{b)y.U{l)^y.U{l)^ 
SM X C/(l)/3. Many studies of Z' are focusing on 
the two extra U{1)' which occur in the above de- 
composition of the i?6. The lightest Z' is defined 
as : 

Z' = Z'^cosp + Z'^sinP (19) 

where the values /? = and /? = 7r/2 corresponds 
to pure Z'-^ and Z,0 states of the x~ ^-iid ■(/'"model, 
respectively. The value /3 = arctan{— ^Jb/i) is re- 
lated to a Z' boson that would originate from the 
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direct breaking of Eq to a rank-5 group in super- 
strings inspired models. 

• the Left-Right Symmetric (LRSM) model is based 
on the symmetry group SUc{3) (8) SUl{2) ® 
SUr{2) (g) U{1)b-l Ifiil, in which B and L are 
the baryon and lepton numbers, respectively. The 
model necessarily incorporates three additional 
gauge bosons and Z' . The most general Z' 
is coupled to a linear combination of right-handed 
and B — L currents: 

J'^R = clb J^j, - (l/2aLfl) J^_^ (20) 



where ulr = \/ {cw9%) / Swdl) - 1> with gL = 
e/sw and gR are the SU{2)l and SU{2)r cou- 
pling constant with = 1 — — siri^Qw- 
The a^^-parameter is restricted to be in the range 
•\/2/3 < ulr ^ \/2. The upper bound corresponds 
to the so-called manifest LRSM with gi — gR, while 
the lower bound corresponds to the x-™odel dis- 
cussed above, since 5*0(10) can results to both 
SU{5) X U{1) and SU{2)r ® SU{2) l «) U{1) break- 
ing parameter. To simplify our study, we will use 
further the following standard assumptions: (i) the 
mixing angles are small; (ii) right-handed CKM 
matrix is identical to the left-handed one, and (iii) 
5fl = 9L- 

• in the sequential model (SSM) the corresponding 
Z' boson has the same couplings to fermions as 
the Z of the SM. The Z' could be considered as 
an excited state of the ordinary Z in models with 
extra dimensions at the weak scale. 

• the Stueckelberg extension of the SM (StSM) is 
based on the gauge group SU{3)c SU{2)l ® 
U{1)y <E) U{1)x m. This extension of the SM in- 
volves a mixing of the U{1)y hypercharge gauge 
field and the U{l)x Stueckelberg gauge field. The 
Stueckelberg gauge field has no couplings to the 
visible sector fields, while it may couple to a hid- 
den sector, and thus the new physical Z' gauge bo- 
son connects with the visible sector only via mix- 
ing with the gauge bosons of the physical sector. 
These mixings, however, must be small because of 
the LEP electroweak constraints and consequently 
the couplings of the Z' boson to the visible matter 
fields are extra weak, leading to a very narrow Z' 
resonance. The width of such a boson could be as 
low as a few MeV or even lower. An exploration of 
the Stueckelberg Z' boson in the Tevatron data was 
recently carried out in [65j. Such a resonance may 
also be detectable via the Drell-Yan process at the 
LHC by an analysis of a dilepton pair arising from 
the Z' decay [6a|- The coupling structure of the 
Stueckelberg Z' gauge boson with visible matter 
fields is suppressed by small mass mixing parame- 
ters thus leading to a very narrow Z' resonance. Be- 



low we will assume that the Z' — Z mixing strength 
is e = 0.06 [H. 

The Z' boson partial decay width into a fermion- 
antifermion pair is given by 

T{Z' ^ fj) = Nc^^^/T^j X 

[{l + 2r^f){gif + {l-Arif){gfj,f] (21) 

where Nc is a color factor {Nc — 3 for quarks and 
Nc = 1 for leptons), g{, g{^ are the left- and right-handed 
couplings of the Z' to the SM fermions, a is the electro- 
magnetic coupling constant, which is a ~ 1/128 at the 
Mz' scale, and -Jrff{— ruf/mz') is assumed to be ^ 1. 
The left-handed couplings of the Z' to the SM neutrinos 

9l="-^ + gl = ^ for Ee and 

LRSM models, respectively, while the right-handed cou- 
plings = in both models. The (g^)^ is restricted to 
lie in the range 0.07(/3 ~ tt/2) < {g{f < 0.45 (/3 ~ 0.4) 
for the Eg model [fiS], and in the range 1/8 < (5{)2 < 3/8 
for the LRSM model. The detail discussions of the StSM 
Z' decay modes and couplings can be found in [66| . 

For the decay rate Z' — 77 in the Z' mass range 
Mz' < 5 TeV one may expect 

T{Z' ^ 77) < 1.4 • Mz' [TeV] MeV, (22) 

as it follows from Eas. (|13ll6l) . 

Due to the possible strong dependence of the couplings 
(/s on the boson mass and the discussion at the end of 
Section II, the inequality (|22|) should be viewed as an 
indication of the ballpark of the possible values of the 
width, rather than the firm limit. The same applies to 
Fig. 3 and caption to Table 1. In Table 1 expected proper- 
ties of Z' bosons for several models are summarized. The 
Z' — > 77 decay rate is calculated by taking into account 
the LEP limit on the Bose violation coupling ([16]). 



model 


Ttot/Mz' 
xlO^ 


B{Z' ^i+fi-) 
xlO^ 


r(Z'^77) 

r(z'^M+M-) 
xlO^ 


ZsSM 


3.1 


3.2 


< 1.5 




0.6 


4.1 


< 5.8 


Zrj 


0.7 


3.5 


< 5.8 




1.4 


5.6 


< 1-7 


Zlrsm 


2.2 


2.2 


< 2.8 


ZstSM 


0.006 


12.3 


< 23 



TABLE I: Summary of expected properties of Z' bosons for 
several models. The first column shows the ratio of the Z' 
width to its mass Mz' , the second column shows the dimuon 
branching fraction, and the third column gives the upper limit 
for the ratio of the Z' diphoton and dimuon decay rates calcu- 
lated taking into account the limit of p6p . For interpretation 
of results, see comments after Eq. (22). 
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IV. SIGNAL AND BACKGROUNDS 

Next, let us explain how to search for the Z' — > 77 
decay at the LHC. The Z' decay into two photons is a 
rare decay mode, with a branching fraction Ji(Z' — )■ 77) 
in the range of < 10^"* — 10^'^. The final state con- 
sists of two high -pT photons, pt — /2 , detected 
in a LHC electromagnetic calorimeter surrounding the 
pp collision region (60 . Isij . The experimental signature 
of the Z' — > 77 decay is a peak in the invariant mass 
distribution of these photons over the continuum back- 
ground. An important point is that if the decays of the Z' 
into leptons occur, the position of the Z' mass peak and 
its expected experimental width in the diphoton invari- 
ant mass distribution can be predicted by the analysis 
of on-peak data based on the observation of a leptonic 
Z' decay mode. The allowed maximal branching fraction 
r(Z' -> 77)/r(Z' -)> At+/Li") calculated for the Bose vi- 
olating coupling constant gsi = 1.0 x 10~^ taking into 
account ([^^ is shown in Table HI One can see, that, 
for example, for the gs — 10~^, the process Z' 77 
could amount to more then 1% of the total muonic Z' 
decay rate in the StSM model. Hence, if Z' is observed 
at the LHC the search for its diphoton decay mode is of 
great interest for possible observation of Bose symmetry 
violation. 




1 1.5 2 2.5 3 3.5 4 4.5 5 

M^. (TeV) 

FIG. 1: The production cross section (j{pp — >■ Z')-B{Z' — >■ 77) 
at the LHC at — 14 TeV in several Z' models indicated in 
the plot , calculated as a function of the Z' mass assuming the 
value of the Bose violating coupling constant gsi = l.Ox 10~^. 

For the appropriate Z' coupling constants discussed 
above, the production cross section a{pp — >■ Z') ■ B{Z' — >• 
77) at the LHC for several Z' models are calculated in the 
framework of PYTHIA 'g^. In this evaluation, the de- 
fault CTEQ5L parton distribution functions [6^] are used 



with no K-factors included. In most of our analysis we 
also neglect errors associated with imprecise knowledge 
of parton distribution functions; the related systematic 
errors, however, will be included into the final results: 
see Sec. VI. Fig[l] shows the production cross section 
cj{pp Z') ■ B{Z' 77) at the LHC at ^ = 14 TeV 
in the Eg, LRSM and Stiickelberg models calculated as 
a function of the Z' mass assuming the value of the cou- 
pling constant gsi — 1.0 x 10^^. The StSM curve is 
calculated for the Z' — Z mixing strength value e — 0.06. 
Interestingly, although the StSM production cross sec- 
tion u{pp — >■ Z') ■ B{Z' — > l^l^) is an order of magni- 
tude below those of other Z' models, the cross section 
cr{pp — > Z') ■ B{Z' — > 77) is comparable to the corre- 
sponding cross sections in other models. 

Although the experimental signature of the decay 
Z' 77 at the LHC is expected to be relatively clean, 
in order to discover this process one has to determine 
whether the Z' mass peak in the diphoton invariant mass 
spectrum could be distinguished from the background 
due to the standard model reactions. At a hadron col- 
lider experiment the diphoton production with a large 
invariant mass is a well known and studied background 
not only for the search of the two photon decay of the 
Higgs boson, but also for searches of new heavy reso- 
nances, extra spatial dimensions, or cascade decays of 
heavy new particles [g^, where it is a source of signif- 
icant background. The dominant standard model back- 
ground sources to our signal are (see e.g. j60l|): 

• the prompt 77 production cither form the quark 
annihilation or gluon fusion. As the final states 
from our signal and from these processes are iden- 
tical, this is irreducible intrinsic background. 

• The 7-l-jets production consisting of two parts: (i) 
prompt photon from hard interaction plus the sec- 
ond photon coming from the outgoing quark due to 
initial and final state radiation and (ii) prompt pho- 
ton from hard interaction plus the decay of a neu- 
tral hadron (mostly isolated) in a jet, which could 
fake a decay photon. The 7 -I- jet, pp 7-f jet + A 
with a jet faking photon production turned out to 
be one of the most important backgrounds. 

• The background from QCD hadronic jets, consist- 
ing of quarks that fragment into a high momentum 
t:^ , which subsequently decays as tt'^ — >■ 77. The 
resulting photon showers may overlap, and can pass 
the photon selection. 

• other possible sources of background is the Drell 
Yan productions pp — > e~'"e~ -I- X. The produc- 
tion of a high-transverse-momentum lepton pair, 
can lead to a diphoton final state if both electrons 
produce hard bremsstrahlung photons, or if the 
electron tracks fail to be properly reconstructed. 
We, however, assume a high efficiency of a LHC 
detector silicon tracker veto which, being applied 
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to photon-containing candidate events, makes this 
background neghgible. 



V. SIMULATIONS OF THE PROCESS 

pp ^ Z' ^ 77 AT ^ = 14 TEV 

To make quantitative estimates, we performed sim- 
phfied simulations at the generator level of the Z' pro- 
duction followed by the decay Z' — 77 in the reaction 
•pp — > Z' and the corresponding background processes 
at the LHC. We consider, as an example, the CMS de- 
tector [6C|. As the signal events preferentially populate 
the large transverse momentum part of the phase space, 
events were generated with Pt > 100 GeV (CKIN(3) pa- 
rameter) and < 2.7 respectively. This allows us to 
reduce the time of computations and also to exclude of a 
very large fraction of the standard model events, which 
are peaked at small transverse momenta. 

The CMS detector is described in detail in Ref. [60| . 
It consists of several subsystems: a superconducting 
magnet, a Si-tracker surrounded by an electromagnetic 
calorimeter followed by a hadronic calorimeter and muon 
chambers used for the detection and reconstruction of the 
events. The CMS experiment uses lead tungstate crystals 
for the electromagnetic calorimeter (ECAL). Each crystal 
measures about 22 x 22 mm^ and covers 0.0175 x 0.0175 
(about 1° ) in the A77 — Acj) space (0 being the azimuth 
angle). 

For photon reconstruction at the generator level, we 
have used the "hybrid" clustering algorithm, to account 
for also fake photons arising from jets [zO] ■ Photon can- 
didates are reconstructed as superclusters in the CMS 
electromagnetic calorimeter (ECAL), within the fiducial 
regions of the barrel (EB) |?7| < 1.4442 and endcaps 
(EE) 1.566 < I77I < 2.5. The superclusters are ex- 
tended in (j) to recover the energy deposited by electron 
bremsstrahlung and photon conversions. We consider a 
photon in the ECAL as a local deposition of electromag- 
netic energy by electrons or photons contained in a cone 
R = y^yy^ + A^ < 0.09 with no associated tracks. 
This definition is equivalent to 10 x 10 crystal size in 
the CMS detector. The CMS experiment uses 5x5 crys- 
tal size to form an energy cluster to reconstruct a photon 
candidate. However, in our efforts to mimic this recon- 
struction process at the generator level, we choose to be 
conservative and use only a 10 x 10 crystal. The momen- 
tum of the photon candidate is defined as the vector sum 
of the momenta of the electromagnetic objects in such a 
crystal. 

As mentioned above, the main challenge to identify- 
ing the true photon candidates arises from jets faking 
photons, see e.g. [G^I- This occurs when a jet from 
the standard model processes with 7 -I- jet or jet — jet 
final state is dominated by a neutral hadron, such as, 
for example, a tt*^ or ry, which decays into two photons. 
If the hadron is highly energetic, so that the cosine of 
the opening angle between the two decay photons is 



cos^Qj^j^) > 0-9, this angle is difficult to resolve and 
the photons can be misidentified as a single energetic 
photon. To suppress such backgrounds, we use various 
isolation variables, without, however, taking into account 
such photon object characteristics as the lateral and lon- 
gitudinal electromagnetic shower shape. Jets typically 
have a larger number of charged particles reconstructed 
in their vicinty, and also a larger ratio of hadronic to elec- 
tromagnetic deposited energy than photons. Likewise, 
hadronic and electromagnetic deposits arising from jets 
will be less isolated than for photons. Fake photon sig- 
nals arising from a jet can be rejected by requiring either 
the absence of charged tracks above a certain minimum 
transverse momentum(Py'^-„ ) associated with the pho- 
ton or the absence of additional energetic particles in an 
annular cone (Riso) around the photon candidate. Fol- 
lowing the diphoton analysis similar to Ref. [7l| , we have 
considered two variables for the isolation purposes: (i) 
the number of tracks (Ntr) from charged particles, such 
as t:, K,p, inside a cone around the photon and (ii) 
the scalar sum of transverse energy (Et) inside a cone 
around the photon. To identify the photons from the de- 
cay Z' — )■ 77, the following 7 events selection criteria are 
used: 

• P^^ > 100 GeV, P}^ > 100 GeV; 

• |77^i^^2| < 2.5, |r;Ti'T2| _^ 1.4442- 1.5666; 

• cos{Qji^j2) < 0.9; 

,Ntr^O for PJiT > 3.0 GeV within Ri^o < 0.35; 

• Et < 5.0 GeV within Riso < 0.35; 

Using this algorithm and requiring the photon to be iso- 
lated, the estimated probability of a jet faking a photon 
in j+jet channel is ~ 10"'' — 10~^. The major sources of 
fake photons are tt" (> 80%), with the rest coming from 
other sources. 

As the next step, the signal event candidates are se- 
lected by requiring that the final state contain at least 
two or more isolated photons and a jet(s). Events are 
studied in which either both photons are in the barrel 
calorimeter, or one photon is in the barrel and the other 
is in the end cup calorimeter. The diphoton invariant 
mass distribution is calculated for two highest Pt pho- 
tons and histogrammed in bins equivalent to the mass 
resolution. The combined acceptance and selection effi- 
ciency for events with Pt > 100 GeV and for M'^ > 1000 
GeV is found to be > 70%. 



VI. RESULTS 

In Fig. [2] the invariant mass M^j distributions in the 
presence of the standard model background are shown 
for events simulated for the LRSM and Ee{Z^) models 
for the Z' with the mass of 1.5 TeV, the Bose violating 
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coupling qbi = 1-0 x 10~^, and the LHC integrated lu- 
minosity L = 150 fb~^. The additional broad peak and 
long tails below the Z' peak are from a combinatorial 
background due to the wrong choice of photons. For the 
invariant 77-mass > 1 TeV, the background under the 
Z' peak drops quickly. 
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from the search for Z — 77 decay mode at LEP, thus 
making the process Z' — > 77 feasible for observation at 
the LHC. For the higher integrated luminosity of 1000 
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FIG. 2: The diphoton invariant mass M-y-y distribution in the 
presence of the SM background calculated for the LRSM and 
'E6{Zx) models for the Z' with the mass of 1.5 TeV and the 
LHC integrated luminosity L = 150 fb'^ . The Bose violating 
coupling is assumed to be gsi ~ 1.0 x 10^^. 



FIG. 3: The estimated CMS discovery potential in the 
{qb, Mz') parameter space calculated for the LRS model and 
for the integrated LHC luminosity 100 and 1000 fb"\ For 
interpretation of results, see comments after Eq. (22). 



The significance of the discovery of the Z' — > 77 events 
can be estimated as [t^ : 

S^2{^/n^Tn^-^/n^), (23) 

where Ns and Nb are the numbers of signal and back- 
ground events respectively, which pass the selection cri- 
teria described above. These numbers of events are esti- 
mated from a search for the Z' — 77 mass peak, which 
was performed in the following way. For every Z' mass 
value, the region around it in the M^^ distribution was 
fitted with a parametrized signal shape centered at the 
My-y value and superimposed over a polynomial back- 
ground. The normalization of each component is allowed 
to float in the fit. This procedure is also used for the 
background estimate as a function of M^^, with statisti- 
cal uncertainties propagated from the fits. The system- 
atic errors discussed below are also propagated to the fit 
procedure. The discovery potential of the Z' -> 77 decay 
with the CMS detector is estimated assuming S >3. 

The final results of this analysis are presented in Fig (31 
where the most stingent limits for the coupling constant 
gB 1 obtained for the LRS model are shown as a function 
of Z' mass. For the total luminosity of 100 fb^^ and 
the mass about 1 TeV, the limits are about one order 
of magnitude lower than the corresponding present limit 



fb ^ , the use of pp Z' reaction allows to probe the 
decay Z' — >■ 77 for the Z' masses up to 2 TeV. 

The systematic errors in the number of expected Z' — > 
77 events coming from various background uncertainties 
are small since the background itself is rather small and 
the discovery region is usually limited by the fast drop 
of the signal cross section at high Z' mass. The largest 
systematic uncertainties of the expected number of Z'- 
bosons arise from the luminosity measurement (6%) and 
the choice of PDF (5%). The latter uncertainty is de- 
termined from the variation in the efficiency when em- 
ploying different PDF parameterizations. To study the 
effect of the detector energy resolution on this analysis, 
the energy of the photons was smeared with the stochas- 
tic term of the CMS electromagnetic calorimeter energy 
resolution ^70.] . Due to limitations in computing time, we 
did not fully simulate the background from jet-jet events. 
Although the dijet cross section is quite large, given the 
low probability of a jet faking a photon it is found that 
the kinematical and isolation cuts used above reject the 
dijet background substantially [7lj. To get conservative 
estimate, we include uncertainties in the background es- 
timate of 15%. Finally, the error of the measured LHC 
integrated luminosity is taken to be 3% [70j . 
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VII. CONCLUSION 

To summarize, we consider a phenomenological model 
of the Bose statistics violation. We show that if a new 
heavy Z' boson is observed at LHC, further searches for 
the Z' — 77 decay mode would suggest an interesting 
additional direction to probe Bose symmetry violation 
at the high energy frontier. We have demonstrated that 
the discovery regions in the (55 1; M^^) parameter space 
for the Z' — >■ 77 decay substantially extend the excluded 
region from the CERN LEP. The low-energy experiment 



in atomic spectroscopy might be a sensitive probe of Bose 
symmetry violation that is complementary to collider ex- 
periments. 
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